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We have fitted low- and medium-energy benchmark datasets employing methods 
used in the MAID/SAID and dynamical model analyses. Independent fits from 
the Mainz, RPI, Yerevan, and Kharkov groups have also been performed over the 
low-energy region. Results for the multipole amplitudes are compared in order to 
gauge the model-dependence of such fits, given identical data and a single method 
for error handling. 



1 Overview 

The following report summarizes results from a series of fits to selected 
datasets for pion photoproduction, a project initiated by the Partial- Wave 
Analysis Working Group of BRAG (the Baryon Resonance Analysis Group). 
The goal of this work was an evaluation of the model-dependence inherent in 
multipole analyses of photoproduction data. In the past, groups have con- 
structed their own databases, and the resulting differences have been shown 
to significantly effect some multipoles, in particular the E2/M1 ratio for the 
A(1232). The handling of systematic errors has also differed, the common 
choices being to either ignore them entirely, to combine them in quadrature 
with statistical errors, or to use them in 'floating' the angular distributions. 

The construction of low-energy (180-450 MeV) and medium-energy (180- 
1200 MeV) datasets was carried out mainly for practical reasons. Many groups 
have studied the region below 450 MeV, which is dominated by the A(1232). 
Only a few groups have extended their analyses over the full resonance region. 
As the recent MAID and dynamical model fits (see the first contribution) 
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extend to 1 GeV, an upper limit of 1.2 GeV was chosen to include as many 
independent fits as possible. p 

The fitted data are listed on the BRAG websitetl and include differential 
cross section, target asymmetry (T), and photon asymmetry (S) data from 
proton targets only. The low- and medium-energy datasets were constructed 
to contain mainly recent measurements, with the goal of minimizing redun- 
dancies and simplifying the fitting procedure. In order to further simplify the 
exercise, systematic errors were not included in the fits. This should be kept 
in mind when values are quoted. 

Below we have compiled the reports of individual groups, giving details of 
the methods used, comments on the fit quality, and ways these results could 
be interpreted. Finally, in the last section, we summarize the findings of this 
study and suggest ways it could be improved or extended. 

2 Multipole Analysis with MAID and a Dynamical Model 
[ S.S. Kamalov, D. Drechsel, L. Tiator, and S.N. Yang ] 

During the last few years we have developed and extended two models for 
the analysis of pion photo and electroproduction, the Dynamical Model □ 
(hereafter called Dubna-Mainz-Taipei (DMT) model) and the Unitary Isobar 
Model H (hereafter called MAID). The final aim of such an analysis is to shed 
more light on the dynamics involved in nucleon resonance excitations and to 
extract N* resonance properties in an unambiguous way. For this purpose 
as a testing ground we will use benchmark data bases recently created and 
distributed among different theoretical groups. 

The crucial point in a study of TV* resonance properties is the separation 
of the total amplitude (in partial channel a = J}) 

.a ,B,a ,.Rm /in 

''777 ''77r I '•777 \^} 

in background t^;^ and resonance t^'^ contributions. In different theoretical 
approaches this procedure is different, and consequently this could lead to dif- 
ferent treatment of the dynamics of N* resonance excitation. As an example 
we will consider the two different models: DMT and MAID. 

In accordance with Ref.H, in the DMT model the t;^^" amplitude is defined 

as 

t^f'{DMT)^e'^'' cosS. 
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where 5a{qE) and R^"^ are the t:N scattering phase shift and full -kN scat- 
tering reaction matrix, in channel a, respectively, qE is the pion on-shell mo- 
mentum. The jpion photoproduction potential v^^"' is constructed in the same 
way as in RefE and contains contributions from the Born terms with an en- 
ergy dependent mixing of pseudovector-pseudoscalar (PV-PS) ttNN coupling 
and t-channel vector meson exchanges. In the DMT model v^^" depends on 

7 parameters: The PV-PS mixing parameter (see Eq.(12) of Ref.H), 4 
coupling constants and 2 cut-off parameters for the vector mesons exchange 
contributions. 

In the extended version of MAID, the s, p, d and / waves of the back- 
ground amplitudes t^^" are defined in accordance with the K-matrix approx- 
imation 

i^^"(MAID) = exp (iS^) cos (Sa<r (<?, W, Q^) , (3) 

where W = E is the total ttN cm. energy and = — fc^ > is the square 
of the virtual photon 4-momentum. Note that in actual calculations, in order 
to take account of inelastic effects, the factor exp (iSa) cos da in Eqs.(2-3) 
is replaced by ^[t^q exp (2i(5q) -I- 1], where both the ttN phase shifts Sa And 
inelasticity parameters rja are taken from the analysis of the SAID groupQ. 

From Eqs. and one finds that the difference between the back- 
ground terms of MAID and of the DMT model is that pion off-shell rescat- 
tering contributions (principal value integral) are not included in the back- 
ground of MAID. From our previous studies of the p-wave multipoles in the 
(3,3) channel □ it follows that they are effectively included in the resonance 
sector leading to the dressing of the jNA vertex. However, in the case of s 
waves the DMT results show that off-shell rescattering contributions are very 
important for the £'o+ multipole in the n'^p channel. In this case they have to 
be taken into account explicitly. Therefore, in the extended version of MAID 
we have introduced a new phcnomcnological term in order to improve the 
description of the vr*^ photoproduction at low energies, 

AE 

E,orr{MAID) = ^^^j^2q2^y Pom , (4) 

where Fd is the standard nucleon dipole form factor, B — 0.71 fm and AE 
is a free parameter which can be fixed by fitting the low energy tt'^ photopro- 
duction data. Thus the background contribution in MAID finally depends on 

8 parameters. Below Tr+n threshold for both models we also, take into account 
the cusp effect due to unitarity, as it was described in RefE, i.e. 

+ 




Ecusp = —o-TrN Rei^J^ W 1 J , (5) 



combined: submitted to World Scientific on February 8, 2008 3 



where cu and cUc = 140 MeV are the tt"*" cm. energies corresponding to W = 
Ep + Ej and Wc = m„ + TO^+, respectively, and a^^^ — 0-12A/m^+ is the pion 
charge exchange amplitude. _ 

For the resonance contributions, following Ref.l3, in both models the Breit- 
Wigner form is assumed, i.e. 

where f-^ji is the usual Breit-Wigner factor describing the decay of a resonance 
R with total width Tb,{W) and physical mass M/j. The phase (t>R{W) in 
Eq. is introduced to adjust the phase of the total multipole to equal the 
corresponding ttN phase shift 5^- 

The main subject of our study in the resonance sector is the determination 
of the strengths of the electromagnetic transitions described by the amphtudes 
A^{Q^). In general, they are considered as free parameters which have to be 
extracted from the analysis of the experimental data. In our two models 
we have included contributions from the 8 most important resonances, listed 
in the Table 1. The total number of amplitudes is 12 and they can be 
expressed also in terms of the 12 standard helicity elements A1/2 and ^3/2- 
Thus, to analyze experimental data we have a total of 19 parameters in DMT 
and 20 parameters in MAID. The final results obtained after the fitting of the 
high-energy (HE) benchmark data base with 3270 data points in the photon 
energy range 180 < E^ < 1200 MeV are given in Table 1. 

For the analysis of the low-energy (LE) data base with 1287 data points 
in the photon energy range 180 < E^ < 450 MeV in the DMT model we 
used only 4 parameters: The PV-PS mixing parameter and 3 parameters 
for the P33(1232) and Pii(1440) resonances. In MAID we have one more 
parameter due to the low energy correction given by Eq. (^. The fi- 
nal results for the helicity elements and the E2/M1 ratio (REM) are given 
in Table 2. In Table 3 we summarize our results and show the oh- 
tained for different channels and different observables after fitting the LE 
and HE data bases. Note that the largest in the LE fit we get for 
differential cross sections and target asymmetries in p(7,7r'^)p. Similar re- 
sults were obtained practically in all other analyses. A detailed compari- 
son with the results of different theoretical groups is given on the website 
http : / / gwdac.phys.g'wu.edu/analysis/prJ)enchrnark.htrn\ Below, in Fig. 



1 we show only one interesting example, the £^0-1- multipole in the channel 

with total isospin 1/2. In this channel contributions from the 5*11(1535) and 

5*11(1620) resonances are very important. At E.y > 750 MeV our values for 
1/2 . 

the real part of the p-Eo+ aniphtude are mostly negative and lower than the 
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N* 




MAID 

current 


MAID 
HE fit 


DMT 
HE fit 


PDG2000 


P33(1232) 


^1/2 


-138 


-143 




-135± 6 




A3/2 


-256 


-264 





-255± 8 


Pi 1 (1440) 


l/Z 


-71 


-81 


-77 


-65± 4 


£>13(1520) 


Al/2 


-17 


-6 


-7 


-24± 9 




A/2 


164 


160 


165 


166± 5 


5ii(1535) 


Al/2 


67 


81 


102 


90± 30 


S3i(1620) 


Al/2 





86 


37 


27± 11 


Sii(1650) 


Al/2 


39 


32 


34 


53± 16 


Pi5(1680) 


Al/2 


-10 


5 


10 


-15± 6 




A3/2 


138 


137 


132 


133 ± 12 


£>33(1700) 


Al/2 


86 


119 


107 


104± 15 




A3/2 


85 


82 


74 


85± 22 


PV-PS mixing: 


Am 


450 


406 


302 




AE 


2.01 


1.73 






xVd.o.f. 


11.5 


6.10 


6.10 





Table 1. Proton hclicity amplitudes (in 10~^ GeV^/^), values of the PV-PS mixing pa- 
rameter Am (in MeV) and low-energy correction parameter AE (in 10~^/m^+) obtained 
after the high-energy (HE) fit 



N* 


MAID 

current 


MAID 
LE fit 


DMT 
LE fit 


PDG2000 


P33(123 2) Al/2 
A3/2 


-138 
-256 


-142 
-265 




-135± 6 
-255± 8 


Pn(1440) A/2 


-71 


-81 


-93 


-65± 4 


REM(%) 


-2.2 


-1.9 


-2.1 


-2.5± 0.5 


xVd-o.f. 


4.76 


4.56 


3.59 





Table 2. Proton helicity elements (in lO'^ GeV'^^'^) and REM=E2/M1 ratio (in %) ob- 
tained from the LE fit. 



results of the SAID multipole analysis. The only possibility to remove such a 
discrepancy in our two models would be to introduce a third Su resonance. 

1 /2 

Another interesting result is related to the imaginary part of the pi?o+ 
plitude and, consequently, to the value of the helicity elements given in Table 
1. Within the DMT model for the 511(1535) we obtain ^1/2 = 102 for a total 
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Observables 


LE 

N MAID DMT 


HE 

N MAID DMT 




317 4.68 3.32 
354 7.22 5.74 


871 6.36 5.95 
859 6.87 5.85 


E(7,^+) 
5](7,7r") 


245 2.79 2.57 
192 2.22 1.58 


546 4.57 6.49 
488 7.65 7.65 


T(7,7r+) 
T(7,7r0) 


107 3.28 2.94 
72 5.18 4.84 


265 3.75 4.17 
241 5.31 5.65 


Total 


1287 4.56 3.64 


3270 6.10 6.10 



Table 3. jN for the cross sections (^), photon (E) and target (T) asymmetries in 
(7, 7r+) and (7, tt") channels obtained after LE and HE fit. N is the number of data points 




150 300 450 600 750 900 1050 1200 150 300 450 600 750 900 1050 1200 
E^(MeV) E^(MeV) 

Figure 1. pE^^^ multipolc obtained after the HE fit using MAID (solid curves) and DMT 
(dashed curves). The dash-dotted curves and data points are the results of the global and 
single-energy fits obtained by the SAID group. 

width of 120 MeV, which is more consistent with the results obtained in 77 
photoproduction, than with previous pion photoproduction results obtained 
by the SAID and MAID groups. 

3 Multipole Analysis of Pion Photoproduction with 

Constraints from Fixed-t Dispersion Relations and Unitcirity 

[ O. Hanstein, D. Drechsel and L. Tiator ] 
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3.1 Outline of the Analysis 

The method presented in Ref. i has been used to analyze the benchmark data 
set. 

The starting point of our analysis is the fixed-t dispersion relation for the 
invariant amplitudea3, 



Re4.(s,t) = 4'P°'<=(s,i) 



1 



V / ds' 



ImAi{s',t), (7) 



with the Mandelstam variables s, u and t, the isospin index / = 0, ±, Sthr = 
(mjv + ni-^Y, = ±1 and = ±1. The pole terms ^^'^"'''(s, t) are obtained 
by evaluating the Born approximation in pseudoscalar coupling. 

The multipole projection of the dispersion relations (0) leads to a system 
of coupled integral equations of the form 



7 

KeM\ {W) = TWf + -V dW ^ /Cf,, {W, W')lmMl (W) 



(8) 



where A4j (W) denotes any of the multipoles El_^_ or M;{|_. The integral kernels 
IC[i,{W,W') are regular kinematical functions except for the diagonal kernels 
/Cf;, which contain a term cx 1/{W — W). 



15 



7 
o 



2. 5 



-1 r 



100 




200 300 400 
(MeV) 



500 




200 300 400 
E, (MeV) 



500 



Figure 2. The real and imaginary parts and the pole term (p.t.) contributions of the 
amplitudes _Eo+|rt/2) and i?o+(3/2). The results of our fit (solid lines) are compared with 
those from Ref. E3 (dashed lines). The data points are the result of our energy independent 
fit. 
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After appropriate simplifications, namely neglect of weakly coupling inte- 
gral kernels and restriction to the partial waves of low angular momentum, the 
integral equations (H) can. be solved in different ways. We based our analysis 
on the method of Omnesu because this method leads to a parametrization of 
the multipolcs in a natural way and thus allows for analyzing experimental 
data. As an additional input, the solution of the integral equations requires 
the phases 0[ (W) of the multipoles on the whole range of integration. Accord- 
ing to the Fermi- Watson theorem, these phases are equal to the corresponding 
ttN scattering phase shifts below 27r threshold, <j>j{W) — d{{W). Above 2tt 
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5.0 
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E, (MeV) 
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7 
o 
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real y lj\ 
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1 
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Figure 5. The real and imaginary parts of the ampHtudcs A/i+(l/2) and Mi+(3/2). Sym- 
bols as in Fig. ^. 



threshold, we use the ansatzes 

(hj (W) = arctan , ' \ ' \ ' . 

for the Sii, Pii, P33, and D13 waves, and 

tjI(W) sm2(5f(W^) 



4'i (^) — arctan 



l + Vi(W)cos26{{W) J ' 



(9) 



(10) 



for the 5*31, P13, and P31 waves. These ansatzes are each based on unitarit 
and an additional assumption. They contain the scattering phases shifts 
Sj (W) and the inelasticity parameters rij {W) ofirN scattering. Both ansatzes 
give (W) — S{{W) below 2tt threshold. Since partial wave analyses of ttN 
scattering are only available up to about W = 2 GeV, we cut the integrals 
off at this energy to avoid the integration of unknown functions. Instead we 
represented the contributions of the imaginary parts at higher energies by 
i-channel exchange of vector mesons. The integral equations then take the 
form 



hl*{W')Mj{W')dW' 
W -W-is 



(11) 



A 
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where A = 2 GeV. 

The solutions relevant for our case are the sum of a particular solution, 
which contains the inhomogeneities as driving terms, and a solution to the 
homogeneous equation multiplied by an arbitrary real coefficients cj : 

Ml{W)^Ml^''^'\W)+cJiMi'''°"\W). (12) 

The coefficients cj are the fitting parameters in our procedure. In addition, 
we varied the coupling constants of the co and the p meson. Since the addition 
of the homogeneous solution is only allowed for multipoles for which the phase 
is different from zero in the asymptotic limit, we decided to fit the parameters 
cf only for the following multipoles: Eo+{0), Eo+il/2), Mi_(0), Mi_(l/2), 
£'i+(3/2) and Mi+(3/2). So we end up with a 10 parameter fit, the results 
of which are discussed in the next section. 



3.2 Results 

The fit to the benchmark data set leads to an overall of 3.7 per data point. 
This high value is mainly due to the differential cross sections and target 
asymmetries of tt" production (see Table |). 



Table 4. The por data point in our fit for the individual observables. It is seen that, 
except for the beam asymmetry S, the description of tt" production data is much poorer 
than that of 7r+ production data. 





7r"p TT'^n 


da/dfl 

S 

T 


7.4 2.8 
1.2 1.6 
5.2 2.9 



The magnitudes of the vector meson coupling constants as determined 
by our fit differ somewhat from the values quoted in the literature. This can 
be attributed to the fact that the dispersion integrals up to 2 GeV already 
contain a certain fraction of the vector meson contributions (see Table p). 

Our results for the s- and p-wave multipoles am, shown in Figs. 3 to 
^. As the procedure presented by Inna Aznauryar£3 at this workshop is 
closely related to ours, her results are also shown for comparison. The real 
parts of the i?o+ and A/i_ multipoles differ significantly from the pole term 
contributions. In our approach, this difference, which in other approaches has 
to be provided by mechanisms like rescattering, is due to contributions from 
dispersion integrals. 
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Table 5. Comparison between the vector meson coupling constants resulting from our fit 
and values quoted in the literature. 





9p gi 9}!> gi 


this 
Rcf. 
Rcf. 


m 
11 

12 


prk 


4.85 15.69 6.78 -1.67 
3.24 19.81 15.85 
1.99 12.42 20.86 -3.41 



4 Analysis of Low-Energy Benchmark Data Using Fixed-t 
Dispersion Relations 

[ I.G.Aznauryan ] 

In this analysis the low-energy benchmark data are analyzed using fixed- 
t dispersion relations Hrit.hin the approach which is close to the approach 
developed in Refs. Ejll3'E3ll3. The real parts of the amplitudes are constructed 
through real parts of invariant amplitudes A^^"^\s,t)^ i — 1,4 (Ref. 0). 
which are obtained using fixed-t dispersion relations: 

ReA^^'Hs^t) = Af^^^{s,t) + ^ 7 W(±^°)(.',t) + ds', 

n J \s' — s s— u J 

Sthr 

(13) 

with 

°'^(s, t) = Af °™(,s, t) + A^s, t) + AP{s, t) (14) 

where -q^^'^^ = —■q'^^^ = 1, ?7i = ?72 = — ?73 = ?74 = 1- Af°''"^{s,t) are the 
contributions of the nuclcon poles in the s- arui. u-channels and of the pion 
pole in the t-channel. In the dispersion relationsll3, the dispersion integrals are 
taken over resonance energy region up to s^ax — (2 GeV)"^ , and it is supposed 
that the integrals over higher energies can be approximated by the t-channel 
u) and p-contributions: A^Ls,t), A^[s,t). These contributions are taken in 
the form presented in Ref. □. The integrals over resonance energy region are 
saturated by the resonances used in the VPI analysis of pion photoproduction 
dataa. The coupling constants for all resonances, except P33(1232), arc taken 
from this analysis. The resonaaee contributions are parametrized in the Breit- 
Wigner form according to RefJla. 

3/2 3/2 

The multipole amplitudes M^j^ and E^j^ corresponding to the Ei^ll2'i2) 
resonance are parametrized using the approach developed in Refs. IjO. Ac- 
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52.458 ±0.002 mFm 


Ini£:i(^ 


-1.19 ± 0.01 mFm 


E2/MI 


-0.023 



3/2 3/2 

Table 6. The amplitudes M^j^ resonance position. 





Present analysis 


RefB 


Ref.l 


Ref.§ 




4 


6.8 


21 


8-14 




-24 


-1.7 


-12 


0-(-14) 




2.5 


4.9 


2 


1.8-3.2 


9e 


21 


16 


13 


8-21 



Table 7. The loNN and qNN coupling constants in comparison with results from other 
sources. 



3/2 3/2 

cording to this approach the amphtudes M-^'j^ , E-^^'j^ are the solutions of the 
singular integral equations which follow from dispersion relations for these 
amplitudes, and have the form: 

M{s) = M^J-{s) + M-,^,{s) + cmM''-(s), (15) 

where M{s) denotes any of the amplitudes M'^^ , ^part^{s) and 

Mp^^((s) are the particular solutions of the integral equations generated by 
the Born and ld contributions. Particular solutions have definite magnitudes 
fixed by these contributions. M^°™{s) is the solution of the homogeneous 
part of the integral equation; it has a certain energy dependence fixed by the 
integral equation and an arbitrary weight which was found by fitting the data. 

In the P33(1232) resonance region we have taken into account also the 
contributions of the non-resonant multipoles E'q^^, £^9+, £-0+^, , M^Jj* 

into lmA^^'^\s,t). These multipole amplitudes were found by calculating 
their real parts from the dispersion relationsll3; then the imaginary parts of 
the multipole amplitudes at < 1.3 GeV were found using the Watson 
theorem. At higher energies the smooth cutoff of these contributions was 
made. 



Our fitting parameters were: 

(1) the constants cm and ce which correspond to the magnitudes of ho- 
mogeneous solutions for M'^^ , in Eq. 15 at the resonance position; 
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cr(7r+) 


A(7r+) 


T(7r+) 


(T(7r") 






X2/N 


922/317 


440/253 


245/107 


2568/347 


276/192 


394/71 



Table 8. The values of x2. 



(2) the coupling constants gY, , gZi 9^ ^ 5 J which describe the lo and p 
contributions in Eq. 14. 

(3) we have also included into the fitting procedure the coupling constants 
for the resonances S'ii(1535), Pii(1430) and Di3(1520), namely, the cou- 
pling constants for the multipoles , e'^^ of 511(1535), for the multipoles 
mI'^, M^°} of Pii(1430), and for the multipoles E^'^ , E^°l, M^L^, M^°} of 
1513(1520). A small variation oLthese coupling constants around the values 
obtained in the analysis of Ref. □ was allowed. 

The obtained results are presented in Tables 6-8. 

5 The EfTective Lagrangian Analysis of the Benchmark Dataset 

[ R. M. Davidson ] 

5.1 The Model 

Details of the effectiveJLagrangian approach(ELA) to pion photoproduction 
may be found in Ref.E2l. Here I briefly summarize the main features of the 
model. The effective Lagrangian consists of the pseudovector (PV) nucleon 
Born terms, i-channel to and p exchange, and s- and u-channel A(1232) ex- 
changes. At the tree-level, the amplitude is gauge invariant, Lorentz invariant, 
crossing symmetric, and satisfies the LET's for these reactions to order m^/M. 
However, the tree-level amplitude violates unitarity. To unitarize this ampli- 
tude, the tree-level multipoles, A4f, are projected out and unitarized via a 
K-matrix approach; 

Ml = Mf cos Sie'^' , (16) 

where 6i is the appropriate nN elastic scattering phase shifts. In practice, 
this is done only for the s- and p-wave multipoles. In order to keep multipoles 
of alL I values, the unitarized multipoles are added to the tree-level CGLN 
J^'sc3 and the tree-level multipoles are subtracted. To give a simple example, 
if only _Bo+ is unitarized, then J-i in this model is 

= + El+ (cos Joe'*" - 1) , (17) 
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where TJ is the tree- level approxiiiiatioii to Ti. 

After unitarization, one is now ready to fix the parameters of the model. 
The A mass, Ma, and width F are determined by a fit to the P33 phase shift 
and are not varied in the photoproduction fit. The pion and nucleon masses 
are fixed at 139.6 and 938.9 MeV, respectively. Although one could look at the 
sensitivity of the photoproduction data to the pion-nucleon coupling constant, 
I keep the PV coupling constant, /, fixed at a value of 1.0. The p and u) are 
taken to be degenerate with a mass of 770 MeV, and the Vir^ {V = p ov u>) 
coupling constants are taken from the known radiative decays V ^ nj. The 
Dirac and Pauli-like VNN coupling constants were allowed to vary in the 
fit. The remaining parameters of the model are related to the A interactions. 
Of most interest are the two 'yNA coupling constants, gi and g2, which are 
related to Ml and E2 by 



where Aja = (M^ — M^)/(2Ma). The remaining parameters are the off-shell 
parameters, X,Y,Z, associated with the A transition vertices. For example, 
the ttNA Lagrangian is of the form 



where a depends linearly on Z. Thus, the off-shell parameters essentially con- 
trol the relative strength of the "ffi'fv term compared to the g^i, term. It 
should be noted that when calculating matrix elements involving the off-shell 
parameters the pole in the A propagator is canceled. Thus, the contributions 
involving the off-shell parameters appear as contact terms. As the off-shell 
parameters arc fitted to the data, these contact terms can partially compen- 
sate for the lack of strong form factors which might arise, for example, from 
pionic dressing of the vertices. 

5.2 Results and Discussion 

The nine parameters of the model were fitted to the low-energy benchmark 
dataset consisting of 1287 data points. The total was 5203 giving a x^/df 
of 4.07. The breakdown of the x^ according to observable is given in Table 
9. It is seen that the highest X^/i^ for da/dQ{Tr°) and T(7r°), which was 
true for the other analyses also. The best agreement was with the photon 
asymmetry, S, for both and 7r+. 
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C - K^{g^u + a-ff,-f„)Nd''TT + h.c. , 



(18) 



As the in this analysis is sUghtly larger than in the other analyses, it 
is useful to look at the breakdown in energy bins. This is shown in Table 
10, where [x, y) and [x, y] have their usual mathematical meaning. It is seen 
that the fit to the lowest energy bin is extremely poor. Most of the data in 
this energy interval are 7r° differential cross section data, and a comparison of 
the fit with these data shows that the fit has too large of a forward-backward 
asymmetry as compared to the data. At these energies, this asymmetry is 
determined by the interference between the Eq^ and the p-wave multipoles. 
A comparison of this multipole obtained in this fit with the same multipole 
obtained from fits that reproduce these low-energy data shows a significant 
difference. 



Table 9. for each observable fitted. 



OBS 


n 






da/dn{-K+) 


317 


988.4 


3.1 


da/dn{Tr°) 


354 


2804.0 


7.9 


S(7r+) 


245 


388.0 


1.6 


E(7r") 


192 


258.1 


1.3 


T(7r+) 


107 


274.9 


2.6 


T(7rO) 


72 


489.6 


6.8 



Table 10. 3s a function of energy bin fitted. 



Interval 


n 






[180,210) 


64 


747.7 


11.7 


[210, 240) 


137 


767.6 


5.6 


[240, 270) 


208 


803.6 


3.9 


[270, 300) 


181 


578.3 


3.2 


[300, 330) 


191 


516.8 


2.7 


[330, 390) 


157 


598.6 


3.8 


[360, 390) 


120 


265.5 


2.2 


[390, 420) 


132 


307.6 


2.3 


[420, 450) 


97 


615.4 


6.3 



At the BRAG workshop, we were able to understand this diff'erence. As 
L. Tiator pointed out, at low-energies for tt" production, the coupling to the 
TT+n channel can be important and one must somehow account for processes 
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like 

7P — > titt"'" — > pvr" . (19) 

In dispersion relation modelsE and dynamical modelsi, this channel coupling 
is dynamically taken into account. In partial wave analysesu, the paramatriza- 
tion is flexible enough to account for this physics. In the isobar modeld, this 
important physics was included in a semi-phenomenological manner. In the 
ELA, this channel coupling is only partially taken into account via the uni- 
tarization. However, dispersive corrections are not explicitly accounted for. 
It is implicitly assumed that the main effect of the dispersive corrections is 
to renormalize the parameters to their physical values. It is further assumed 
(or hoped) that any additional dispersive corrections can be accounted for by 
the contact terms coming from the off-shell parameters. Evidently, this is not 
the case over the entire fitted energy range. Thus, to improve the low-energy 
fit, without destroying the fit near the peak of the resonance, this additional 
physics would need to be added by hand to the ELA. 

The parameters from the. fit are shown in Table 11 along with represen- 
tative parameters from Ref.Eil. The resulting resonance couplings are 

Ml = 286.2 X lO^^GeV"^/^ 
E2 = -7.21 X lO^^GeV"^/^ 

^ = -2.55% . (20) 
Ml ^ ' 



Table 11. Parameters obtained in this fit (B.M.) compared with those obtained in 



Par. 


B.M. 


DMW 




5.06 


5.01 ± 0.22 


.92 


4.72 


5.71 ± 0.43 


Z 


-0.25 


-0.30 ± 0.12 


Y 


1.65 


-0.38 ± 0.66 


X 


-5.74 


1.94 ± 2.28 


.9^1 


-0.09 


8 




4.46 


-8 


9pi 


1.51 


2.66 


gp2 


3.66 


16.2 



The gi coupling, which is mostly responsible for the Ml strength, has 
not changed much compared to the earlier work. g2, which determines the 
strength of E2 has changed within the error bar. As a rule of thumb, the 
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smaller 32, the larger in magnitude E2. The claaijige in 32, and hence E2, is 
probably due to the new high-precision E data c3c3. The ofF-shell parameter 
Z has not changed much, but the other parameters have. This is partly due to 
the fact that the vector meson couplings were allowed to vary in this fit, but 
not in RefO. Since 1 did not do an error analysis for the BM fit, it is hard to 
say if the changes are significant. I do know that the vector meson couplings 
and off-shell parameters are highly correlated fit parameters. I should point 
out that the vector mesons are put in as point particles, i.e., no form factors. 
Thus, when I fit the vector meson coupling 'constants', I would expect them 
to decrease compared to their on-shell values. However, the results for the ui 
seem unusual. The general conclusion at the BRAG workshop is that the A 
region is not a good place to fit the vector meson couplings. 

From a glance through the various multipole solutions, it seems to me 
that the model dependence is under good control in the A region. One mul- 
tipole that remains to be understood is the M]^_ . Though all analyses pretty 
much agree on its numerical value, its physical interpretation is quite differ- 
ent in the various models. Is it a crossed-A effect or the tail of the Roper? 
In principle, the amplitude should be crossing symmetric, so if ther|e_is an s- 
channel A exchange, there must be a u-channel A exchange. In Ref.cil, it was 

1 /2 

found that the Mj^i multipole is largely insensitive to the off-shell parame- 
ters. Thus, these contributions cannot suppress the u-channel contribution in 
this multipole, and there is no room for a Roper contribution_as large as in 
the isobar model. What is even more mysterious is that in Ref.Eil it was found 
that the Roper contributiO|LL.enters with the opposite sign than the u-channel 
contribution (see Fig. 8 in Ell). Regarding the Roper contribution to this mul- 
tipole, there are two clear possibilities. Either someone has made a mistake, 
or two different things are being compared. Certainly the latter is true to 
some extent. In Ref.cil, both the s- and the u-channel Roper exchanges were 
included, whereas in the isobar modeS only the s-channel contribution was 
included. Thus, it is possible that the s- and u-channel contributions destruc- 
tively interfere in this multipole to produce the small effect found in c3. I do 

not know the solution to this puzzle, but if dispersion relations are telling us 

3/2 

there is a large crossing contribution to this multipole from the A/j^^ , then we 
must necessarily investigate this problem within the framework of a crossing 
symmetric model. 
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6 Multipole Analysis of the Benchmark Set Covering the First 
Resonance Region 

[ A.S. Omelaenko ] 

6. 1 Introduction 

A great deal of experimental data on the photoproduction of single pions is 
stored in compilations. There are problems connected with the normaliza- 
tion of some systematic measurements of the differential cross section. As a 
result, in multipole analyses, one has to reject or renormalize up to 10% of 
experimental points. This can throw some doubt on the extraction of delicate 
values, such as the E2/M1 ratio for the A"*" (1232). This makes very interest- 
ing the comparison of analyses, based on different approaches, using the same 
test database for the most well investigated 7p — > p7r''(n7r+) reaction, with a 
strong accent on the modern data. Here we report on a fit to the low-energy 
dataset based on (a) the resonance model with polynomial parametrization of 
the background and (b) energy-independent fitting. 

6.2 Formalism in the First Resonance Region 

Resonance Model. Description of the P33 nN scattering amplitude is sim- 
ple if it is taken to be purely elastic in the first resonance region. Similar 
to the form used in Ref.L3, our multi-parameter model for single pion photo- 
production is written as a sum of a background plus the A+(1232) resonance 
contribution. The real part of the background is given by the electric Born 
approximation which is completed in s-, p- and d- waves by a cubic polynomial 
of the form 

4 4 
ReM4(£;,) = ReM4(i?W) J] {E, - E^^^)/{E^^ - (21) 
1=1 3=1 

Here, depending on the laboratory nhpton energy E^^ multipoles Mfj^ = 
, are defined according to Ref.EJ with the isospin structure 

A^/"^ = l/iA^Op + \/2/3A^+„, 

Index I is the orbital angular momentum, while -f and — correspond to the 
total angular momentum j ^ I ± 1/2. Imaginary parts of the background 
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multipoles are calculated according to Watson's theorem, using phase shifts 
from TriV elastic scattering analyses: 

ImM/i = ReM/i tan(52/,2(;±) ) • (22) 

3 /2 3 /2 

The resonant multipole amplitudes A^_^ and B^_^ have the foUowing form 

(^^ ) = Mf (S^ ) + Sm {E^ ) cos <533e'''^^ (23) 
with the first term describing excitation of the resonance 



^1+ (^7) - ^Af;/^::r7^72 T772 TF^' (24) 



with W and Wq being the total cm. energy and its value at resonance, 
respectively. Cm being the resonance constant. The energy-dependent widths 
were parametrized by 



r{w)^ro{^] (26) 



where k and q are the cm. momenta of the photon and pion, respectively, 
and fco, qo are the corresponding values at Wq. The second tepn in Eq. 
(23) corresponds to Noelle's unitary treatment of the backgroundnj. In our 
approach, the resonance term has a simplified form, not containing the elastic 
background phase shift 6b- However, this has practically no influence on the 
E2/M1 ratio, being quite independent of (5bE3- The real background functions 
Ba{E^) and Bb{E^) are parametrized according to Eq. (21), each in terms 
of 4 knot values. The phase shift ^33 is calculated as the phase of the Breit- 
Wigner form in Eq. (24): 

''''''' -W^' ^ ^ 

The resonance quantities CA^ C'b, Wq, Tq were determined in the fit, along 
with the knot values of Ba and Bb and the real parts of the background 
multipoles: Ao+, Ai+, A2-, B2- with / = 1/2 and Ao+, Ai_ 

with / = 3/2. Inclusion of the d-wave multipoles is motivated by a possible 
influence of the second resonance region. 

Energy-Independent Version. In order to fit narrow energy bins, we 
take the real parts of the abovementioned non-resonant multipoles as indepen- 
dent parameters. Eq. (22) is used to calculate imaginary parts, taking into 
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account the energy dependence of the phase shifts. To avoid calculational 
problems at the point where 633 passes through tt/2, the resonant multipoles 
have been parametrized as 

=^exp(*J33), (28) 
B^{'f'> ^ BcM^Sss), (29) 

with A and B to be determined for each energy bin along with the real parts 
of all background multipoles. 



6.3 Results from the Fits 

Resonance Model. The whole energy interval was split into 3 subintervals, 
with knot values corresponding to the central values for energy-independent 
fits. 



^^,2,3,4)^200,250,300,350 MeV 
^a,2,3,4) ^250,300,350,400 MeV 
E:^^'^'^'''' =300,350,400,450 MeV 



Fit 1 
Fit 2 
Fit 3 



180< <350 MeV 
250< <400 MeV 
300< E^ <450 MeV 



Fit 


Data 


a;" 


X 


aVn 


Fit 1 


Total 


2026.5 


930 


2.18 




7P — » pir^ 


1332.4 


460 


2.90 




7P — > n7r+ 


694.1 


470 


1.48 


Fit 2 


Total 


1759.2 


871 


2.02 




7P pir^ 


1115.7 


388 


2.88 




7P — > riTr"*" 


643.5 


483 


1.33 


Fit 3 


Total 


1259.7 


697 


1.81 




7P — > pn^ 


821.6 


306 


2.69 




7P n7r+ 


438.1 


391 


1.12 



Table 12. Fitting in the resonance model framework. 



For each fit a total of 44 free parameters were searched (the parameter 
X turned out to be difficult to determine, and was fixed at Walker's value 
of 185 MeV). The statistical characteristics for all fits are given in Table 12, 
separately for the neutral and charged particle production. In Table 13, the 
contribution to from individual experiments is given, as calculated from 
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T^p'r\ ft ion 


Observable 


Label 


N 


A / 


IP ^ 




FTIQfiMA 


45 


2.2 




rlrr MO 




^9 
oz 








BE97MA 


77 


1.4 




Sigma 


BL92LE 


16 


3.3 




Sigma 


BE97MA 


77 


0.8 




T 


B098BO 


28 


3.1 




T 


BE83KH 


26 


3.6 


— ^ n7r+ 


da/dfl 


FI72BO 


32 


0.9 




da/dQ 


BROOMA 


39 


1.9 




da/dfl 


BEOOMA 


140 


1.7 




Sigma 


BLOOLE 


57 


1.7 




Sigma 


BEOOMA 


140 


0.9 




T 


DU96BO 


75 


1.1 



Table 13. Fits with the resonance model in the 250-400 MeV interval. See the database 
summary webpage to match abbreviations with full references. 



Interval (MeV) 


Mass (MeV) 


Width (MeV) 


E2/M1 (%) 


180-350 


1232.6 ± 1.0 


121.1 ± 6.0 


-1.83±0.31 


250-400 


1231.7± 0.7 


118.4 ± 3.9 


-2.34±0.13 


300-450 


1230.6 ± 0.8 


116.5 ± 3.9 


-2.26±0.14 



Table 14. Resonance parameters for the A(1232) from different energy intervals. 



Fit 2. Finally, in Table 14, values for the resonance mass and width are given 
along with the corresponding value for the E2/M1 ratio. 

Energy-Independent Analysis. Results from the resonance model fits 
were used as starting values for 10 parameters, which were determined by 
minimization from the data contained in each 10 MeV bin. For all formally 
successful fits, the corresponding central values for photon energy, number of 
points (N) and per degree of freedom is given in Table 15. 

6.4 Concluding Remarks 

In the low-energy region, dominated by the first resonance, we find a rather 
good determination of the main s- and p-wave partial-wave amplitudes for 
pion photoproduction (on proton targets), taking into account d-wave correc- 
tions. Differences between the Fits 1-3 on overlapping energy intervals can be 
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il'-y , iVlC V 


X 


1\T 
IN 


X /"/ 


Xl/-y , iVlC V 


X 


IN 


X /"J 


210 


10.1 


18 


1.27 


330 


100.6 


51 


2.45 


220 


18.6 


43 


0.56 


340 


20.0 


49 


0.51 


230 


23.8 


22 


1.99 


350 


118.3 


74 


1.85 


240 


52.9 


58 


1.10 


360 


6.3 


32 


0.29 


260 


112.3 


71 


1.84 


370 


28.9 


42 


0.90 


270 


33.3 


36 


1.28 


380 


93.1 


46 


2.59 


280 


171.5 


69 


2.91 


390 


16.7 


42 


0.52 


290 


72.9 


51 


1.78 


400 


86.8 


56 


1.89 


300 


81.3 


80 


1.16 


410 


16.7 


34 


0.70 


310 


121.8 


59 


2.49 


420 


67.2 


50 


1.68 


320 


96.8 


69 


1.64 











Table 15. Statistics of the energy-independent fits. 



considered as some measure of the model error. In our analysis, a suspicious 

1/2 1/2 

deviation in the energy dependence of the M■^^_ and M^^^ was evident over 
the 380-420 MeV interval. Encouraging was the rather stable determination 
of the mass and width of the A"*" (1232) resonance and E2/M1 ratio within 
the framework of this model (Fits 2 and 3). 

A comment of the chosen database is in order. Of the cross section data 
for neutral pion production, the set of HA97MA gives a which is definitely 
too large. This is true in all of the fits. Problems seem concentrated mainly 
at the small and large angles. Unfortunately, the modern polarization data 
do not represent the complete set of single polarization observables. As a 
result, in some bins of the energy-independent fit, the experimental data were 
not sufficient for a successful minimization procedure, or yielded unreasonably 
low (less than 1) values for x^ /df . Further examination has shown that the 
abovementioned deviations in the energy behavior of the d-wave multipoles 
are due to the lack of reliable data as well. Precise new experiments would be 
extremely desirable. 

7 SAID Multipole Analysis of the Benchmark Dataset 
[ R.A. Arndt, I.I. Strakovsky, R.L. Workman ] 

SAID fits were performed on both the low- and medium-energy datasets, using 
a phenomenological method described in Ref.Q. Multipoles were parametrized 
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in the form 

M = (^phen + Born) (1 + iT^n) + B^^^^T^n (30) 

where Aphcn and -Bphcn were constructed from polynomials in energy having 
the correct threshold behavior. An additional overall phase proportional to 
ImTjrjv — ^^AT ^^'^ allowed in order to allow deviations from the above 
form for waves strongly coupled to channels other than tvN . This form clearly 
satisfies Watson's theorem, where valid, and allows a smooth departure as 
new channels become important. A real Born contribution is assumed for 
unsearched high partial-waves. 

The overall fit to data is summarized in the Table below. While the 
fit to most quantities is about 2 per data point, there are clear exceptions. 
Our fit to the neutral-pion differential cross section and target polarization is 
significantly worse in the low-energy region. If is calculated for the region 
between 450-1200 MeV a fairly uniform fit emerges, with all data types fitted 
equally well. 



Observable 


Low-Energy Fit 


Medium-Energy Fit 


Medium-Energy Fit 




(180-450)MeV 


(180-1200)MeV 


(450-1200) MeV 




4.69 


3.14 


1.80 




1.24 


2.05 


2.32 




4.49 


2.84 


1.89 




2.12 


2.55 


2.08 




1.46 


1.81 


2.36 




2.92 


2.25 


1.83 



Table 16. Comparison of the low- and medium-energy fits in terms of x^. 



The poor fit to neutral-pion difi'erential cross sections can be traced to 
the Haerter data from Mainz (Ph.D. Thesis, unpublished - see the database 
description on the BRAG webpage). This set contributed a x^/datum of 10, 
the remaining data again having an overall per data point of about 2. 

Why the fit to target polarization should be poor is harder to determine. 
In the low-energy region, the benchmark dataset has taken target-polarization 
data from Kharkov and Bonn. Of these, the fit to Kharkov data is particularly 
bad. Our fit to the Bonn data is reasonable apart from the lowest-energy (and 
perhaps also the highest energy) set. 

8 Conclusions and Future Projects 
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8.1 Database Issues 



As this exercise was motivated by an earlier study of the E2/M1 ratio, and 
its database dependence, we should (at least) expect to verify the relative 
model- independence of this quantity (Issshen a standard dataset is selected). 
This point was discussed by DavidsonEa in his contribution to NSTAR2001. 
From the combined set of benchmark fits, the E2/M1 ratio was found to be 
—2.38 ± 0.27%. Here the important quantity is the error, since a change in 
the fitted dataset could shift the central value. 

Qualitative features of the fit to data were reasonably similar in the low- 
energy region. There was general agreement that the differential cross section 
and target polarization were badly fitted in the neutral-pion channel. In 
all cases, the thesis data of Haerter were identified as problematic for the 
differential cmss section. [These data can be abandoned, as a new Mainz 
measuremento of the cross section and photon asymmetry for neutral-pion 
photoproduction has been analyzed and is nearing publication. This new set 
covers the full angular range and a wide range of energies.] 

Unfortunately, such a simple solution was not evident for the target po- 
larization. The fit to both the Bonn and Kharkov data is poor, with some 
sets suggesting a different shape. In this case, the general disagreement should 
motivate a re-measurement of this quantity, hopefully with better precision. 
A confirmation of the shape suggested by these data would be very difficult 
to accommodate. 

8.2 Model Dependence 

One ingredient common to all analyses is the Born contribution, usually aug- 
mented with vector-meson exchange diagrams. While it is tempting to deter- 
mine an optimal set of vector meson couplings, over the fitted energy range, 
our study has shown this to be a very unstable procedure, when fits are re- 
stricted to the low-energy database. 

1/2 

The M^_ multipole, and its low-energy behavior, provides an interesting 
example of the model dependence encountered when one tries to decompose 
a partial-wave into its underlying ingredients. Hanstein has noted that, in his 

3/2 1/2 1/2 1/2 

dispersion integral, the coupling M^^^ exceeds M^__ . In 

the study of Davidson, this is accounted for by a u-channel A exchange, a 
contribution absent in many isobar-model approaches. As a result, the Roper 
resonance contribution depends strongly on the set of approximations defining 
a fit. 
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8.3 Further Comparisons 

One unique feature of this study was the construction of a sitd3 containing, in 
one place, both the database and the multipole fits. This facihty was buih into 
the SAID siteQ, thus aUowing a wide range of comparisons. Users can compare 
fits to observables (in both the benchmark and fuU database), the resulting 
multipoles, and search for regions of maximal agreement or disagreement. 



8.4 Extensions 0/ this Study 

Having the present set of fits as a guide, it would be useful to repeat this study 
over a more carefully chosen database - containing, in some way, the effects 
of systematic errors. A second project of interest would be the extension to 
electroproduction, with the goal of understanding the model dependence seen 
in the extracted Ei_|_/Mi+ and S1+/M1+ ratios as a function of . 
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